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We work out the dynamics of the compressible edge of the quantum Hall system based on the 
electrostatic model of Chklovskii et al. We introduce a generalized version of Wen's hydrodynamic 
quantization approach to the dynamics of sharp edge and rederive Aleiner and Glazman's earlier 
result of multiple density modes. Bosonic operators of density excitations are used to construct 
fermions at the interface of the compressible and incompressible region. We also analyze the dy- 
namics starting with the second-quantized Hamiltonian in the lowest Landau level and work out the 
time development of density operators. Contrary to the hydrodynamic results, the density modes are 
strongly coupled. We argue that the coupling suppresses the propagation of all acoustic modes, and 
that the excitations with large wavevectors are subject to decay due to coupling to the dissipative 
acoustic modes. A possible correction to the tunneling density of states is discussed. 



PACS numbers: 73.40.Hm, 73.20.Mf 



I. INTRODUCTION 



The theory of edge excitations in the fractional quantum hall liquid had been first derived by WenS based on the 
gauge invariance of the Chern-Simons term for a bounded system. A more intuitive picture of the edge excitation was 
^ laid out in his subsequent paperl where a classical energy and a classical equation of motion for the incompressible 
liquid is written down, and he stipulated that the conjugate variables in the corresponding Lagrangian satisfy the 
canonical commutation relations. In the absence of a more microscopic theory, such "phenomenological quantization" 



,: serves as a guideline in our understanding of the phenomena and not infrequently yields accurate results. Qwite a lot 
' of progress in the study of vortex motion in superfluid He and superconductors had been made in this wayjj'cl 
\C Eventually, any phenomenological theory has to be based on the analysis of the underlying microscopic Hamiltonian. 

A good example is the plasma oscillation in a metal where one can show that a quantum- mechanical eauation of motion 
of the density operator reduces to the (semi-)classical one in the regime where RPA is applicable.l3 The strength of 
the microscopic analysis lies in the fact that it yields more information than can be garnered at the phenomenological 
Ch ' level. Coupling of a plasma mode to single-particle excitations give rise to Landau damping and yields a correction 

. to classical dispersion that can be systematically calculated using many-body theory. 
^ ' With a sharp edge, there are so few degrees of freedom that one cannot distiH^aush between collective modes 
Q , and single-particle excitations. This leads to bosonization of one-dimensional fluidiloEl for which the single-particle 
O properties can be uniquely derived once we know the collective excitation spectrum. 

^ The situation changes when one realizes that the physical edge is often not a sharp boundary separating an 

• incompressible bulk from its surroundings. Chklovskii et alii have demonstrated that with the gate-confined 2DEG, 
the edge is a rather wide, smoothly varying region of density. While the idea of compressible edge had been actively 
adopted both experimentally and theoreticallya in the past few years, no attempt, to our knowledge, has been made 
to bridge it with the bosonization theory of the edge. A number of worksphave sprung up recently to address the 
effect of long-range Coulomb interaction on the properties of bosonized edgeja but the geometry is still that of a sharp 
edge as in Wen's original paper|,-pne interesting attempt to understand the dynamics of the compressible region was 
made by Aleiner and GlazmanE2l, henceforth referred to as AG, who used classical hydrodynamics of charged fluid 
to derive the (necessarily collective) spectrum of the compressible edge analytically. They showed there are many 
density modes within the compressible edge, one of which is charged and higher in energy than the other, so-called 
acoustic modes, which are neutral. Their model, however, does not account for the boundary effect taking place at the 
interface of compressible and incompressible regions. We have tried in this paper to understand the gapless spectrum 
of a wide edge more rigorously from the quantum- mechanical standpoint. 

In Sec. II we will show how the phenomenological quantization approach of Wen can be generalized to the dynamics 
of compressible edge including its boundary with the bulk, and recover AG's result along the way. We also discuss 
the issue of bosonization in this kind of edge. In Sec. Ill, we start from the second-quantized Hamiltonian to address 
the same problem. It turns out that the density modes are strongly coupled because of the diffuse distribution of the 
electron occupation. In Sec. IV, we discuss various consequences of such coupling, including its possible effect on the 
tunneling density of states. 
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II. PHENOMENOLOGICAL APPROACH TO EDGE DYNAMICS 



The geometry of the sample is that of a semi-infinite Hall bar where the boundary separating the compressible and 
incompressible liquid is defined as y = and the compressible strip extends in the positive y direction (fig. 1(a)). 

In the presence of a relatively weak local electric field fluctuation and a strong external magnetic field, one can 
write Vx — SEy/B and Vy = —SEx/B where the field is produced either by external potential or as a result of electron- 
electron interaction. In turn, the local field can be expressed as the gradient of the change in potential energy. To 
completely specify the motion, we need the following equations for the local density n(x, y). 



dn dn dn 
dt dx ^ dy 

n{x,y + h{x,y,t),t) = no{y). (2.2) 



:«^ + — ^y = 0' (2-1) 



The first equation is the continuity equation for the divergence-free velocity field, dvx/dx + dvy/dy = 0. The sec- 
ond one parameterizes the dynamics by identifying points on the density contour no{y) with (x,y + h{x,y)). The 
situation is depicted in fig. 1(b). The boundary between the compressible and incompressible region is defined by 
n{x, h{x,0,t),t) — rio(O) = v/^ttIq where Iq is the magnetic length. We will set both this and Planck's constant, ?i, 
equal to 1 in the rest of the section. If one imagined a collection of very long strings stacked along y-direction, each 
carrying a label no(y), the density fluctuation can be understood in terms of the strings wiggling about equilibrium, 
but never so violently enough that they cross other strings. 

One can deduce a set of identities by differentiating Eq. (2.2) with respect to a;,?/, and i, which can be put back 
into Eq.(pTl) to give 

' ■ ■ ^^TJZ 1 + 7^ =0- (2-3) 



dy \dt dx J \ dy 

We have dno/dy <0 inside the compressible region. For small fluctuations, dh/dy should be small so that it is only 
the middle term which vanishes identically. By defining h(x,y,t) = dri{x,y,t)/dx, one arrives at the equation of 
motion of the compressible edge 

d^rj d'^r] 1 d dri{x,y,t) 



where V{x, y + h) is the local potential(fluctuation) at a given time t. One can also check that varying the following 
action reproduces the above equation of motion. 

' ^ dxdydt — [ — / eV{x,y + y')dy'dxdydt (2-5) 



2 dy dxdt J dy 

with e being the electric charge. Quantization comes from the commutator 



[77(x, y), ^h{x\ y')] = i6{x - x')6{y - y'). (2.6) 



Note that results derived in eqs. ( p.4| )-(2.6) reduce to known expressions for sharp edge by taking the limit dno{y)/dy — 
-{v/2ir)5{y). .... 

When one is solely interested in the dynamics inside the compressible edge, as AG were, one can regard 
V{x,y) as Coulomb potential due to the fluctuating density. In linear approximation, n{x,y) ~ nQ^y) —5n{x,y) ~ 
— (dno/dy) h{x,y), and the Hamiltonian is given, up to quadratic order in density, by 

/■ dnodno h{x,y)h{x' ,y') ,,,/,/ . 
H = — ===dxdydx dy ; k = Aireoe (2.7 

2k J dy dy' ^{x ~ x')^ + {y - y')^ 

One can expand the operator h{x, y) in terms of its Fourier components in xfi <x <L. Each component hq{y) satisfies 

2p2fl 1"°° dvn 

ih,{y) = ^ dy'^K,{\qy - qy'\)h,{y'). (2.8) 

where Kq is the modified Bessel function of zeroth order. For an eigenmode hqa{y) whose time dependence is e~*"°*, 
it becomes a homogeneous integral equation. We can exploit the monotone nature of the function no(y) to use tiq as 
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a coordinate, and remove the partial derivative inside the integral. Hence we arrive at a Fredholm integral equation 
with-a real, symmetric kernel, for which the eigenfunctions form a complete, orthogonal set and all eigenvalues are 
real.llil The kernel Kq is expressed in terms of eigenvalues and eigenfunctions by 

oo 

Koilqyino) - qy{n'„)\) = ^ ^^.(g) V("o) VK)- (2-9) 



a=0 



The eqs. (|2.8[)-(p.9[) summarize the results initially obtained by AG. When we have the boundary separating the 
incompressible and compressible(IC) region, density fluctuation in the compressible region necessarily implies some 
of the electrons on the incompressible side have to move as well. When they move, they are doing work against the 
field that confines them, and it costs additional energy. The amount of energy necessary is given by 

Hic='^ I h{x,O fdx. (2.10) 
47r J 

where y — indicates the IC boundary. Indeed, this was the only term that appeared in the Hamiltonian for a sharp 
edge, without the electron-electron interaction. 

The total Hamiltonian, Htot = H + Hjc is written in the basis of eigenfunctions as 



h{x,y) ^ L ^^^y^ Cqahqa{y)e~ 



iqx 



E ( ^ + -£^2^(0)) ct„c,„ + E ^ V(0)/.,,(0)ct„c,,. (2.11) 



The boundary term contributes a diagonal term which modifies the slope of the dispersion but, more importantly, 
there are off-diagonal terms which mix states that were thought to be orthogonal. The perturbation is quadratic, 
hence by rotating to a proper basis, one can again recover independent modes. One might as well think of Eq. 



(2.9) as containing the complete dynamical in form ation, provided we use the eigenmodes and eigenvectors of the 



full Hamiltonian. Quantization condition, Eq.(2.6), becomes [cJ„,Cgc[] — q and the Hamiltonian becomes a sum of 
oscillator modes whose energy levels are set by ujaiq). Before mixing, ujQ{q) cx q|logq| is the most dominant mode 
while a > 1 are nearly linear. The boundary will mix these and the linear modes become weakly logarithmic. 

FromJjhe set of harmonic oscillator-type operators at our disposal, we can follow the standard recipe of bosonization 
in Idtlflla and construct the following set of fermion operators at the IC boundary. 

i^iix) (X exp{E J^kcq.e-^''' " cl^e^'l} (2-12) 

Now ijjl^{x) involves a coherent density change over the whole compressible region including the IC boundary and 
differs from Wen's construction in which an electron is created at the perimeter of the incompressible fluid. 

III. MICROSCOPIC APPROACH TO EDGE DYNAMICS 

Our starting point is the Coulomb Hamiltonian, 

p2 r 1 

H=— -^UrUUr'U(r'U(r)(frd^r'. (3.1) 

2k J \r ~ r'\ 

We are interested in the dynamics within the lowest Landau level, therefore we can expand tp{r) as follows;@ 

V>(r) = iV^L)-^/^ E e'^"e-5(?/-fc)'cfc. (3.2) 



When it is substituted into Eq. (3.1), we obtain 
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kik2q 



2 

qa 



Pli<l) = — 1=ff XI 'i«"(^)4+g/2Cfc-9/2 = Pa{-q)- (3.3) 



The difference between the normal-ordered form, Eq. (3.1), and the above Hamiltonian is an infinitely large constant 



which can be cancelled by assuming a compensating positive charge background. The Hamiltonian is approximate in 
the sense that we have replaced the convolution f f{y)e~'^y~^^ ^"^dy with \/27r/(fc) or, in other words, identified y with 



Zgfc. We have used the expansion of Kq, Eq. (2_^), and it should be understood that the effect of the IC boundary 



can be included by resorting to appropriate generalizations of u>a {q) and hqa ■ 

More can be understood about the form of the above Hamiltonian by recalling that for a homogeneous system, the 
Coulomb Hamiltonian takes the form i/'2J2^^{Q)p{'l)p{^Q): provided one discards the exchange term. Because the 
y-direction is not translationally invariant, we have a taking over the role of a momentum, labeling the underlying 
classical transverse modes. 

The hydrodynamic result of the previous section is tantamount to having each density operator, Pj^(q), oscillating 
at a definite frequency, iVaiq)- We will examine the Heisenberg equation of motion of the density operators as follows: 

-^H(p) = ^E^"(«){^i(9)[^"(«)'^/^(^')] + [puQ),pMpim- (3.4) 

qa 

We have absorbed the overlap integral e^'^^/^ into uJa{q)- The r.h.s. of the above equation contains the hydrodynamic 
term, ti;^(p)Pg(p), but there are also off-diagonal components present. Evaluation of those terms must rely on some 
kind of approximation scheme. We will solve Eq. (^^) by linearizing the r.h.s. so that it reduces to a set of coupled 
linear equations for Pg(p). 

The first method we try is a version of RPA, in ''jiilich one only keeps those bilinear terms with the momentum 
transfer equal to p, i.e. terms proportional to c^^pCfc.l3't2l To leading order in momenta, we get 

2 

P<t(q)[P„(q),Pj(p)]« 

The primes indicate derivatives with respect to fc. In deriving the above result, we have used 

[Pa((j), Flip)] « (1/ \/\pq\L) ^{ph'q^hpp + qhgah'pp} '^k+{p~q)/2^k-(p-q)/2- (3.6) 

fc 

The diagonal term in Eq. ( |3.5[ ), when multiplied by L0a{q)/2 and summed over q and a, gives half of uj0{p)Pp{p). 
The other half comes from [PQ,(g), Pg(p)]Pj((7). In evaluating the second term in Eq. (3^), we assume that various 



eigenfunctions do not depend on q as much as they do on a so that we can simply label them ha{k). The integration 
over q extends from to cx), but due to rapidly decaying e""^ ' , there is a cutoff at large momenta. Since our interest 



lies in estimation, we conveniently choose the region of integration for q to be [0,p]- The Eq. (3.4) becomes 



^fi{p)Pl{p) -Eu;a(p) ^ ^{hlhf})' nkcl^p/^Ck-p/2 w/3(p)P^(p) - ^ Wa(p)^a/37^'-| (p) (3-7) 

a ^ ' k a, 7 

where Aafi-y is a constant of order Ig/aL, a being the width of the compressible edge. The size of coupling of Pg(p) to 
P^{p) is the sum u;a{p)Aapj ^ Hl/aL) ^a{p) ~ {hlL){'^a{p))avg- In arriving at the last relation, we assumed 
the net number of modes involved are roughly a/^o- So the coupling between modes is generally weak. 

The underlying assumption of the above analysis, or indeed in the original RPA theory, is that there exists a stable 
collective excitation. For a stable excitation, p is a SQod quantum number, and it may be safe to keep just those terms 
that preserve p. In the theory of Nozieres and PineaJ, it is a posteriori justified by calculating various decay rates of a 
plasmon which turns out to be small compared to the energy scale of a plasmon. In our problem, there is no obvious 
energy scale to set collective modes apart from single particle excitations and there may be a large coupling between 
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the two. This may in turn show up as a couphng between modes of different momenta and if it is significant, will 
destroy the assumption of the RPA. 



As another way to linearize Eq. (3.4), we propose to replace the commutators with their expectation values in the 
ground state. The resulting equation will range over differing momenta as well as a. Taking the expectation value of 
Eq. (3J) requires the knowledge of (c^_|_pCfe). For p = 0, it is just the density, Uk- For py^O, we can argue the following. 



If the state k was unoccupied in the ground state, Cfc acting on the ground state would destroy it. Likewise, if the 
state k+p had been occupied, there would be no place for an added electron to go, hence the expectation value will 
be again zero. So a plausible answer is (c|_|_pCfc) = nk{l—nk+p) ~ nk{l—nk) for small p. We get precisely this result 

for the form of the ground state \g.s.) = J^^(\/l— rife + y^nfeC^)|uac). There are as yet no definitive answers in the 
literature on how the ground state of the edge should behave at T = 0. Our guess is that if one is at a sufficiently high 
temperature that any many-body correlations which may exist at absolute zero are destroyed, but still low enough 
that quantum hall effect survives, the ground state should be approximately described by the single particle state 
wavefunction given above. 

The commutator part can now be worked out straightforwardly. For p = q, it is equal to Saf3- When p ^ q, 

([P^iq),P;ip)]) = P^S^, + -f^A^p, (3.8) 

where Sa(3 — Jq hah/3{l — 2nk)dnk, and Aap — {hf^dha / dnk—hadhj3 / dnk){nk—nl)dnk ■ The two integrals contribute 
numbers of order unity as shown in the Appendix. We find it more convenient to carry out the analysis in terms of 
PI (p) I y/p and will refer to it as P^ (p) . We write cua (p)-^ ^aP where it is understood that the velocities have weak 
logarithmic dependence and uo/wi ~ — log(pa) w log lO'^JlJ With v = l and assuming the existence of just two modes, 
a = 0, 1, we have 

- zPj(p) = c^o(p)Pj(p) + ^u;,ip) J2 pHq) 

TT ^ — ' 

-^Plip) = i^,ip)Pl{p) + ^Y.''oiq)P^iq). (3.9) 

9 

What we have now is a coupled equation without the presence of small dimensionless parameters, lo/a, or a/L. This 
is the crucial feature that arises as a result of the microscopic calculation. The modes which were thought to be 
independent in the hydrodynamic model are in fact strongly coupled. One can do a perturbative analysis of the above 
equation using the fact that ll!i{p)/luo{p) — vi/vq is small, 

9 9 

-^P^{p)^cJ,{p)P^{p)+^^\ a;i(p)5]^Pj(g'). (3.10) 

V / 9 9' 

Up to this order of vi/vq, the dynamical behavior of Pl{p) is completely determined by that of Pp (p) and results in 
renormalization of the equation of motion for P^{p). With v = 1/3, there is also coupling within the same mode, 

- zPj(p) = u;,{p)Pl{p) + J2 ^voP^iQ) +J2^'^^(P-1 + '^P/^')Plil) 

9 9 

~^P}{p) = u;^{p)P}{p) + ^v,P}{q) + J2 ^Mq -P + ^Q/^')P^iQ)- (3-11) 

9 9 



Why is there so much mixing? For a sharp edge, the occupation number makes a step function, and we get 



{cl., Ck) = SpoUk- Compressible edge has a diffuse distribution of electrons, giving a fairly large contribution to 



{cl^pCk) for p different than zero, and hence to the coupling of modes. In the light of the above analysis, our initial 
RPA assumption to keep only the momentum-conserving terms seems unjustifiable. 
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IV. DISCUSSION 



The density excitation at the edge is often called the edge ma g netopl asmo n (EMP) and to distinguish the different 



modes, the a > 1 are called the "acoustic" modes. The Eqs. (3^) and ( 3.11 ) then represent the coupling of the EMP 



modes with the acoustic modes. Since experiments are ofte n do ne at the f undamental frequency corresponding to 



p — ^-k/L, one can keep just Pq{2t:/L) and Pi{2ti/L) in Eqs. (3.9) and (3.11) and disregard higher modes. Note that 
the Gaussian factor implicit in the frequency uja tends to suppress the coupling to higher momentum modes. We also 
introduce damping of the acoustic mode, wi(p) = (wi + i"l)p where 7 is comparable or larger than vi but still quite 
small compared to vq. This is the problem of two coupled harmonic oscillators one of which is heavily damped. The 
resulting frequencies are easily found to be VQ^Vi) + \^{vi+i"f) and ^ (1 — A^)(fi + ry), with A = A\/2i' / 1:"^ . As a 
result of coupling, the initially dissipationless mode becomes weakly dissipative, and the acoustic mode propagates 
more slowly. The acoustic modes, a > 1, were predicted to be essentially unobservable due to dissipationEEl ia. 
currently available samples with mobilities 100 m^/Vsec. A number of EMP experiments both in the integertil 
and fractionalllJilj have been unable to detect the diversity of modes either. Our analysis shows the mode coupling 
works in favor of slowing down the speed of propagation of the acoustic mode, thus making its detection more difficult. 
For highe r momentum states, the effect of coupling should be more significant, as one can see from the second half 



of Eq. (3.1C). Since wi is complex, the more terms added in the summation will increase the dissipation. We therefore 
expect the states with large p to decay easily compared to the modes close to the fundamental frequency. 

In the experimental study of Luttinger liquid behavior of the edge, it is necessary that the edge be as sharp as 
possible in order to both reduce tha number of active modes in the edge, and to eliminate coupling between them. 
Recent experiment of Chang et aZJlJ was designed to create a very sharp edge profile, and the experiment seems to 
vindicate many predictions of the chiral Luttinger liquid theory. They see exponents in the I-V measurement which 
are not totally integral, as expected from the theory, but rather show deviations of the order of 10%. In the light 
of this result, it is instructive to examine the possible correction of the tunneling density of states away from the 
ideal sharp-edge, short-range interaction limit. The long-range character of the dispersion makes the evaluati|Qii of the 
spectral function difficult, but it has been successfully solved in the recent work of Ziiliche and MacDonaldEj. They 
conclude the long-range interaction results in the logarithmic correction to the tunneling density of states. We focus 
instead on the consequence of the coupling of modes. 

We assume a much simplified model with a single- mode edge excitation, such as given by Pl^q) whose excitation is 
undamped for p ranging from to Pc- The modes with p>Pc are damped with a constant coefficient, oj =p+ijp{vo — 1). 
The source of damping will be the coupling to the invisible mode as discussed above. The density of states is 
extracted from the imaginary part of the retarded electron Green's function, G fi{xt) = —i{ip{xt)'il;^ {0)) , and the total 
tunneling depsity of states at energy u is proportional to the imaginary part of Gr{x = 0,w). A straightforward 
calculatiorO shows that 

G«(0,c.)(x^'"-W eM^t + m -e'^*'^"^ (1 - g-^^*)} (4.1) 

Jo Jo 9 [t + t-ft ~ tO+)"' 

where m — l/v. For pc ~ (constant rate of damping for all momenta) the tunneling density of states is proportional 
to w™"^, as is the case for a non-dissipative system given by 7 = O. With pc ^ O, the limit of integration explicitly 
depends on uj and the naive scaling result will be modified. For pc/^^ 1, one can check that the first correction term 
goes like (pc + aj)™~^. One should note, however, that other than this kind of dissipative correction and the effect 
of long-range interaction on the spectral density, the compressible edge will not pose an essential modification of the 
original sharp-edge theory. 
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APPENDIX: EVALUATION OF Scp AND A,,p 

One needs a specific expression of hqa in order to evaluate Saf3 and Aaf3. We follow AG and use hgoluk) — 
sJI-Kjv^ hqa{n]S) — 2yJ TT / V cos{T:ant. / v) for a > 1. One can check this is an exact set of solutions once the charged 
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mode, a = 0, is so large in energy that it decouples from all of the lower modes. There is also mixing due to the IC 
boundary as mentioned in Sec. II, and the true eigenstates are a complicated mixture of the above set. We have not 
considered the mixing of states, but do not expect that it will significantly change the result that follows. 
Using the above eigenfunctions, we obtain 

r a = /3, 27r(l-i/) 

Sap = Spa.= { a = 0^/3, {AV2i^/7T/3^)il - (Al) 

and 

( a ^0^(3, 2V27r{(2;.//327r2)(l-(-l)/5)-(l-i.)(-l)/3} 
= -A^fs ={ a ^13^0, 47r(l - v){{a' + /32)/(a2 _ ^2)|(_^)a+/3^ (A2) 
[ (M7r){(/3 - a)/{(3 + af + (/3 + a)/{f3 - af}{l - (-1)"+'^) 

Some specific results for v ^ 1 and v = l/i are 

5oo = = 5ii, 501=8^2/^, v4oi=8v^/7r {v = \) 

5oo = 47r/3 = S*!!, S'oi = 8\/2/37r, Aqi = 4\/27r/3 + 8\/2/37r [v^l/i). 



^ X. G. Wen, Phys. Rev. B 41, 12838 (1990); X. G. Wen, Phys. Rev. B 43, 11025 (1991) 
^X. G. Wen, Int. J. Mod. Phys. B 6, 1711 (1992). 

^ F. D. M. Haldane and Yong-Shi Wu, Phys. Rev. Lett. 55, 2887 (1985). 
* Ping Ao and D. J. Thouless, Phys. Rev. Lett. 70, 2158 (1993). 

^ D. Pines and P. Nozieres, The Theory of Quantum Liquids (Benjamin, Reading, Mass., 1966). 
® M. Stone, Ann. Phys. 207, 38 (1991). 

^D. B. Chklovskii, B. I. Shklovskii, and L. I. Glazman, Phys. Rev. B 46, 4026 (1992). 

*P. L. McEuen et al, Phys. Rev. B 45, 11419 (1992); C. W. J. Beenakker, Phys. Rev. Lett. 64, 216 (1990); A. M. Chang, 
Solid State Commun. 74, 871 (1990); L. Brey, J. J. Palacios, and C. Tejedor, Phys. Rev. B 47, 13884 (1993); S. W. Hwang, 
D. C. Tsui, and M. Shayegan, Phys. Rev. B 48, 8161 (1993); B. Y. Gelfand and B. L Halperin, Phys. Rev. B 49, 1862 
(1994); M. M. Fogler, E. L Levin and B. I. Shklovskii, Phys. Rev. B 49, 13767 (1994); M. R. Geller and G. Vignale, Phys. 
Rev. B 50, 11714 (1994). 

^ Y. Greg and A. M. Finkel'stein, Phys. Rev. Lett. 74, 3668 (1995); K. Moon and S. M. Girvin, preprint (1995). 
1" I. L. Aleiner and L. I. Glazman, Phys. Rev. Lett. 72, 2935 (1994) 

G. Arfken, Mathematical methods for physicists (Academic Press, New York, 1985), p. 890-893. 

C. de C. Chamon and X. G. Wen, Phys. Rev. B 49, 8227 (1994). 

G. D. Mahan, Many-particle Physics (Plenum Press, New York, 1981). 

V. L Talyanskii et al, Phys. Rev. B 50, 1582 (1994); V. I. Talyanskii et al, J. Phys.: Condens. Matter 7, L435-L443 (1995). 
M. Wassermeier et al, Phys. Rev. B 41, 10287 (1990). 
1^ R. C. Ashoori et al, Phys. Rev. B 45, 3894 (1992). 
" A. M. Chang, L. N. Pfeiffer, and K. W. West, preprint. 



U. Ziiliche and A. H. MacDonald, :ond-mat/9607075 



FIG. 1. Density profile given by n of the edge in equilibrium, (a), and out of equilibrium, (b). The shaded region indicates 
incompressible bulk of filling factor v and a is the width of the edge. The constant density no{y) is given by a straight, horizontal 
line in equilibrium, but is displaced by h{x,y) at each point {x,y) in (b). 
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